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Abstract. The physics of charm has become one of the best laboratories exposing the limitations of the naive constituent 
quark model and also giving hints into a more mature description of meson spectroscopy, beyond the simple quark-antiquark 
configurations. In this talk we review some recent studies of multiquark components in the charm sector and discuss in 
particular exotic and non-exotic four-quark systems. 
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More than thirty years after the so-called November revolution [1], heavy hadron spectroscopy remains a challenge. 
The formerly comfortable world of heavy mesons is shaken by new results [2]. This started in 2003 with the discovery 
of the D*q(2317) and Ds\{2A6Q) mesons in the open-charm sector These positive-parity states have masses lighter 
than expected from quark models, and also smaller widths. Out of the many proposed explanations, the unquenching 
of the naive quark model has been successful [3]. When a [qq) pair occurs in a P-wave but can couple to hadron pairs 
in 5- wave, the latter configuration distorts the {qq) picture. Therefore, the 0+ and 1+ [cs] states predicted above the 
DK{D*K) thresholds couple to the continuum. This mixes meson-meson components in the wave function, an idea 
advocated long ago to explain the spectrum and properties of light-scalar mesons [4]. 

This possibility of {csnh) (n stands for a light quark) components in D* has open the discussion about the presence 
of compact (ccnn) four-quark states in the charmonium spectroscopy. Some states recently found in the hidden-charm 
sector may fit in the simple quark-model description as (cc) pairs (e.g., X(3940), 7(3940), and Z(3940) as radially 
excited XcQ^ Xci^ and Xci), but others appear to be more elusive, in particular X(3872), Z(4430)+, and 7(4260). 
The debate on the nature of these states is open, with special emphasis on the X(3872). Since it was first reported 
by Belle in 2003 [5], it has gradually become the flagship of the new armada of states whose properties make 
their identification as traditional (cc) states unlikely. An average mass of 3871.2 ±0.5 MeV and a narrow width 
of less than 2.3 MeV have been reported for the X(3872). Note the vicinity of this state to the D^D*^ threshold, 
M(D°D*°) = 3871.2 ± 1.2 MeV. With respect to the X(3872) quantum numbers, although some caution is still 
required until better statistic is obtained [6], an isoscalar y^'' = 1++ state seems to be the best candidate to describe 
the properties of the X(3872). 

Another hot sector, at least for theorists, includes the {ccnil) states, which are manifestly exotic with charm 2 and 
baryon number 0. Should they lie below the threshold for dissociation into two ordinary hadrons, they would be narrow 
and show up clearly in the experimental spectrum. There are already estimates of the production rates indicating they 
could be produced and detected at present (and future) experimental facilities [7]. The stability of such (QQqq) states 
has been discussed since the early 80s [8], and there is a consensus that stability is reached when the mass ratio 
M{Q)/m{q) becomes large enough. See, e.g., [9] for Refs. This effect is also found in QCD sum rules [10]. This 
improved binding when M/m increases is due to the same mechanism by which the hydrogen molecule {p,p,e^ ,e^) 
is much more bound than the positronium molecule (e+,e+,e^,e^). What matters is not the Coulomb character of the 
potential, but its property to remain identical when the masses change. In quark physics, this property is named^avoMr 
independence. It is reasonably well satisfied, with departures mainly due to spin-dependent corrections. 

The question is whether stability is already possible for (ccnii) or requires heavier quarks. In Ref. [9], a marginal 
binding was found for a specific potential for which earlier studies found no binding. This illustrates how difficult are 
such four-body calculations. 

Besides trying to unravel the possible existence of bound (ccnn) and (ccnn) states one should aspire to understand 
whether it is possible to differentiate between compact and molecular states. A molecular state may be understood 
as a four-quark state containing a single physical two-meson component, i.e., a unique singlet-singlet component in 
the colour wave function with well-defined spin and isospin quantum numbers. One could expect these states not 



being deeply bound and therefore having a size of the order of the two-meson system, i.e., 1. Opposite to that, 
a compact state may be characterized by its involved structure on the colour space, its wave function containing 
different singlet- singlet components with non negligible probabihties. One would expect such states would be smaller 
than typical two-meson systems, i.e., Ar <l. Let us notice that while Aj? > 1 but finite would correspond to a meson- 
meson molecule A^ oo would represent an unbound threshold. Thus, dealing with four-quark states an important 
question is whether we are in front of a colorless meson-meson molecule or a compact state, i.e., a system with two- 
body colored components. While the first structure would be natural in the naive quark model, the second one would 
open a new area on the hadron spectroscopy. 

There are three different ways of coupUng two quarks and two antiquarks into a colorless state: 



[{qiqimqA)] = {|3i2334), I612634)} = {|33):.M66)12} 

[{qmXqiqA)] = {II13I24), I813824)} = {|ll)c, |88)J (2) 
[{qmXqiqs)] = {|li4l23>, I814823)} = {|l'l')c, |8'8%} , (3) 

being the three of them orthonormal basis. Each coupling scheme allows to define a color basis where the four-quark 
problem can be solved. The first basis, Eq. (1), being the most suitable one to deal with the Pauli principle is made 
entirely of vectors containing hidden-color components. The other two, Eqs. (2) and (3), are hybrid basis containing 
singlet-singlet (physical) and octet-octet (hidden-color) vectors. 

To evaluate the probability of physical channels (singlet-singlet color states) one needs to expand any hidden-color 
vector of the four-quark state color basis in terms of singlet- singlet color vectors. Given a general four-quark state this 
requires to mix terms from two different coupUngs, 2 and 3. In [1 1] the two Hermitian operators that are well-defined 
projectors on the two physical singlet- singlet color states were derived, 
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^in,, = (^e+e^) 2(i-|,(n|i,i-),p) ' 

where P, Q, P, and Q are the projectors over the basis vectors (2) and (3), 



Q = I 88),, (88 I , (5) 



and 



Q = |8'8'),c(8'8'| . (6) 

By using them and the formalism of [11], the four-quark nature (unbound, molecular or compact) can be explored. 
Such a formalism can be applied to any four-quark state, however, it becomes much simpler when distinguishable 
quarks are present. This would be, for example, the case of the [nQnQ) system, where the Pauli principle does not 
apply. In this system the bases (2) and (3) are distinguishable due to the flavor part, they correspond to [{nc){cn)\ and 
[{nh){cc)], and therefore they are orthogonal. This makes that the probability of a physical channel can be evaluated 
in the usual way for orthogonal basis [12]. The non-orthogonal bases formalism is required for those cases where the 
Pauli Principle applies either for the quarks or the antiquarks pairs. Relevant expressions can be found in [11]. We 
show in Table 1 some examples of results obtained for heavy-light tetraquarks. One can see how independently of 
their binding energy, all of them present a sizable octet-octet component when the wave function is expressed in the 2 
coupUng. Let us first of all concentrate on the two unbound states, Af > 0, one with 5 = and one with S=l, given 
in Table 1 . The octet-octet component of basis (2) can be expanded in terms of the vectors of basis (3) as explained in 
the previous section. Then, the probabilities are concentrated into a single physical channel, MM or MM* {MM stands 
for two identical pseudoscalar D (B) mesons and MM* for a pseudoscalar D (B) meson together with its corresponding 
vector excitation, D* (B*)]. In other words, the octet-octet component of the basis (2) or (3) is a consequence of 
having identical quarks and antiquarks. Thus, four-quark unbound states are represented by two isolated mesons. This 
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FIGURE 1. Pmm as a function of A£. 



conclusion is strengthened when studying the root mean square radii, leading to a picture where the two quarks and 
the two antiquarks are far away, (x^)'/^ ^ 1 fm and {y^Y^"^ ^ 1 fm, whereas the quark-antiquark pairs are located at 
a typical distance for a meson, {z^Y^'^ < 1 fm. Let us now turn to the bound states shown in Table 1, < 0, one in 
the charm sector and two in the bottom one. In contrast to the results obtained for unbound states, when the octet-octet 
component of basis (2) is expanded in terms of the vectors of basis (3), one obtains a picture where the probabilities 
in all allowed physical channels are relevant. It is clear that the bound state must be generated by an interaction that 
it is not present in the asymptotic channel, sequestering probability from a single singlet-singlet color vector from the 
interaction between color octets. Such systems are clear examples of compact four-quark states, in other words, they 
cannot be expressed in terms of a single physical channel. 

We have studied the dependence of the probability of a physical channel on the binding energy. For this purpose we 
have considered the simplest system from the numerical point of view, the {S, I) — (0,1) cchh state. Unfortunately, this 
state is unbound for any reasonable set of parameters. Therefore, we bind it by multiplying the interaction between the 
light quarks by a fudge factor. Such a modification does not affect the two-meson threshold while it decreases the mass 
of the four-quark state. The results are illustrated in Figure 1 , showing how in the Ag — > limit, the four-quark wave 
function is almost a pure single physical channel. Close to this limit one would find what could be defined as molecular 
states. When the probability concentrates into a single physical channel (PmiM2 1) the system gets larger than two 
isolated mesons [11]. One can identify the subsystems responsible for increasing the size of the four-quark state. 
Quark-quark ((x^)'/^) and antiquark-antiquark ((y^)'/^) distances grow rapidly while the quark-antiquark distance 
{{^Y^^) remains almost constant. This reinforces our previous result, pointing to the appearance of two-meson-like 
structures whenever the binding energy goes to zero. 

In another recent investigation, the four-body Schrodinger equation has been solved accurately using the hyperspher- 
ical harmonic (HH) formalism [12], with two standard quark models containing a linear confinement supplemented 
by a Fermi-Breit one-gluon exchange interaction (BCN), and also boson exchanges between the light quarks (CQC). 
The model parameters were tuned in the meson and baryon spectra. The results are given in Table 2, indicating the 
quantum numbers of the state studied, the maximum value of the grand angular momentum used in the HH expan- 
sion, /iTmax, and the energy difference between the mass of the four-quark state, E^q, and that of the lowest two-meson 
threshold calculated with the same potential model, A^. For the [cchh) system we have also calculated the radius of 
the four-quark state, 7?4^, and its ratio to the sum of the radii of the lowest two-meson threshold, A^. 

As can be seen in Table 2 (left), in the case of the [ccnh) there appear no compact bound states for any set of 
quantum numbers, including the suggested assignment for the X(3872). Independently of the quai^k-quark interaction 



TABLE 1. Heavy-light four-quark state properties for selected 
quantum numbers. All states have positive parity and total or- 
bital angular momentum L = 0. Energies are given in MeV. The 
notation M\M2 \t stands for mesons M\ and M2 with a rela- 
tive orbital angular momentum (.. /'[|33)'^(|66)'^)] stands for 
the probability of the 33(66) components given in Equation (1) 
and P[ \ ll)c (1 88),)] for the 11(88) components given in Equa- 
tion (2). Pmm, Pmm' , and Pm*m* have been calculated following 
the formaUsm of [11], and they represent the probability of find- 
ing two-pseudoscalar {Pmm), a pseudoscalar and a vector {Pmm*) 
or two vector {Pm'W ) mesons 



(5,/) 

Flavor 


(0,1) 
ccnn 


(1,1) 

ccnn 


(1,0) 
ccnn 


(1,0) 
bbnn 


(0,0) 
bbnn 


Energy 
Threshold 


3877 

00 Is 
+5 


3952 
DD* \s 
-1-15 


3861 
DD* Is 
-76 


10395 
BB* \s 
-217 


10948 
BiB\p 
-153 


P[133)12] 

pmYc"] 


0.333 
0.667 


0.333 
0.667 


0.881 
0.119 


0.974 
0.026 


0.981 
0.019 


p[mc] 


0.556 
0.444 


0.556 
0.444 


0.374 
0.626 


0.342 
0.658 


0.340 
0.660 


Pmm 
Pmm* 
Pm'm* 


1.000 
0.000 


1.000 
0.000 


0.505 
0.495 


0.531 
0.469 


0.254 
0.746 



TABLE 2. [ccnn) and {ccnn) results. 
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0++ (24) 
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4441 


-1-15 


0.624 


> 1 


0+- (22) 


4224 


-1-64 


3778 


-1-140 




1+ (24) 


3861 


-76 


0.367 


0.808 


1++ (20) 


3786 


-1-41 


3808 


-1-153 


7 = 


2+ (30) 


4526 


4-27 


0.987 


> 1 


1+- (22) 


3728 


-^45 


3319 


-1-86 




0- (21) 


3996 


-1-59 


0.739 


> 1 


2++ (26) 


3774 


-^29 


3897 


-1-23 




1- (21) 


3938 


-1-66 


0.726 


> 1 


2+- (28) 


4214 


-1-54 


4328 


-1-32 




2- (21) 


4052 


-1-50 


0.817 


> 1 


1-+ (19) 


3829 


-1-84 


3331 


-1-157 




0+ (28) 


3905 


-1-50 


0.817 


> 1 


1— (19) 


3969 


-1-97 


3732 


-1-94 




1+ (24) 


3972 


-1-33 


0.752 


> 1 


0"+ (17) 


3839 


-1-94 


3760 


-1-105 


1 


2+ (30) 


4025 


-1-22 


0.879 


> 1 


0— (17) 


3791 


-H108 


3405 


-1-172 




0- (21) 


4004 


-1-67 


0.814 


> 1 


2"+ (21) 


3820 


-1-75 


3929 


-1-55 




1- (21) 


4427 


-1-1 


0.516 


0.876 


2— (21) 


4054 


-1-52 


4092 


-1-52 




2- (21) 


4461 


-38 


0.465 


0.766 



and the quantum numbers considered, the system evolves to a well separated two-meson state. This is clearly seen in 
the energy, approaching the threshold made of two free mesons, and also in the probabilities of the different colour 
components of the wave function and in the radius [12]. Thus, in any manner one can claim for the existence of a 
compact bound state for the [ccnn) system. 

A completely different behaviour is observed in Table 2 (right). Here, there are some particular quantum numbers 
where the energy is quickly stabihzed below the theoretical threshold. Of particular interest is the l"*" cchh state, 
whose existence was predicted more than twenty years ago [13]. There is a remarkable agreement on the existence of 
an isoscalar 7^ = 1+ cchfi bound state using both BCN and CQC models, if not in its properties. For the CQC model 
the predicted binding energy is large, — 76 MeV, A^; < 1, and a very involved structure of its wave function (the DD* 
component of its wave function only accounts for the 50% of the total probability) what would fit into compact state. 
Opposite to that, the BCN model predicts a rather small binding, —7 MeV, and Ar is larger than 1, although finite. 
This state would naturally correspond to a meson-meson molecule. 



Concerning the other two states that are below threshold in Table 2 (right) a more careful analysis is required. Two- 
meson thresholds must be determined assuming quantum number conservation within exactly the same scheme used in 
the four-quark calculation. Deahng with strongly interacting particles, the two-meson states should have well defined 
total angular momentum, parity, and a properly symmetrized wave function if two identical mesons are considered 
(coupled scheme). When noncentral forces are not taken into account, orbital angular momentum and total spin are 
also good quantum numbers (uncoupled scheme). We would hke to emphasize that although we use central forces in 
our calculation the coupled scheme is the relevant one for observations, since a small non-central component in the 
potential is enough to produce a sizeable effect on the width of a state. These state are below the thresholds given by 
the uncoupled scheme but above the ones given within the coupled scheme what discard these quantum numbers as 
promising candidates for being observed experimentally. 

Binding increases for largerM/m, but in the {bbhh) sector, there is no proliferation of bound states. We have studied 
all ground states of [bbhh) using the same interacting potentials as in the double-charm case. Only four bound states 
have been found, with quantum numbers 7^(7) = l''~(0) , 0"'"(0), 3~(1), and 1~(0). The first three ones correspond to 
compact states. 

To conclude, let us stress again the important difference between the two physical systems which have been 
considered. While for the (ccnn), there are two allowed physical decay channels, (cc) + (nh) and (cn) + (en), for 
the {cchh) only one physical system contains the possible final states, {ch) + [ch). Therefore, a {ccnh) four-quark state 
will hardly present bound states, because the system will reorder itself to become the Ughtest two-meson state, either 
(cc) + {nh) or (cn) + (cn). In other words, if the attraction is provided by the interaction between particles / and j, it 
does also contribute to the asymptotic two-meson state. This does not happen for the {cchh) if the interaction between, 
for example, the two quarks is strongly attractive. In this case there is no asymptotic two-meson state including such 
attraction, and therefore the system might bind. 

Once all possible {cchh), {bbhh) and {ccnh) quantum numbers have been exhausted very few alternatives remain. 
If additional bound four-quark states or higher configuration are experimentally found, then other mechanisms should 
be at work, for instance based on diquarks [4, 14, 15]. 
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